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Abstract. In this paper we develop elements of the global calculus of Fourier 
integral operators in W 1 under minimal decay assumptions on phases and ampli- 
£> \ tudes. We also establish global weighted Sobolev 1? estimates for a class of Fourier 

integral operators that appears in the analysis of global smoothing problems for 
dispersive partial differential equations. As an application, we exhibit a new type 
of smoothing estimates for hyperbolic equations, where the decay of data in space 
is quantitatively translated into the time decay of solutions. 
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1. Introduction 



< 

The paper will deal with a class of Fourier integral operators on MJ 1 that appears 
in the analysis of smoothing problems for dispersive evolution equations, such as 
hyperbolic, Schrodinger, KdV, and others. 

Apart from this motivation, we are interested in global properties of the calcu- 
lus of pseudo-differential and Fourier integral operators under minimal assumptions. 
This follows a general line of research for deriving regularity properties of pseudo- 
differential operators under minimal assumptions on the symbols. Global L 2 bound- 
edness properties of pseudo-differential operators have been thoroughly investigated 
in the literature (e.g. Calderon and Vaillancourt [CVj . Childs |Chj . Coifman and 
Meyer [CM], Cordes |Coj . Sugimoto [Su], etc.). These results have been later ex- 
tended to classes of Fourier integral operators by the authors in |RSj . Local IP 
and other estimates for non-degenerate Fourier integral operators are also quite well 
understood, see e.g. the survey paper by the first author [Rj. 

Results on the global L 2 boundedness of Fourier integral operators in W 1 appeared 
in the papers of Asada and Fujiwara [AFj . Kumano-go |Kuj . Boulkhemair |Boj . with 
different assumptions on the regularity of phases and amplitudes. In all these papers 
the assumption was made that the second order derivatives of the phase with respect 
to the frequency variables must be globally bounded. Unfortunately, this assumption 
fails in applications to global smoothing problems, and it was recently removed by 
the authors in |RS4] . 

The method of dealing with dispersive equations by globally reducing them to nor- 
mal forms has been developed by the authors in |RS2] . This method proved to be very 
effective for a variety of equations. For example, a critical case of the Kato-Yajima 
|KYj smoothing estimate for Schrodinger equations was proved in [RS3J by globally 
reducing the estimate to a model case with useful additional commutator proper- 
ties. The same strategy was later developed in [RS6J to treat smoothing estimates 
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for general dispersive equations with constant and time-dependent coefficients. In- 
stead of previously used spectral methods or dual restriction theorems from harmonic 
analysis, the approach of [RS2j-[RS5j relies on global L 2 estimates for the canoni- 
cal transforms that in turn can be realised as a class of globally defined Fourier 
integral operators in IR n . Moreover, smoothing estimates themselves are essentially 
global weighted Sobolev L 2 estimates for solutions to dispersive equations, and to- 
gether with Strichartz estimates they constitute the main tools for dealing with time 
well-posedness questions for nonlinear evolution equations. 

Transforming smoothing estimates by the method of canonical transforms leads to 
the necessity to analyse weighted Sobolev estimates for Fourier integral operators. 
In |RS4j . the authors investigated questions of L 2 and weighted L 2 boundedness for 
the appearing operators. In Theorem 13.51 of this paper we will now give a result in 
weighted Sobolev spaces over L 2 . This result will rely on a global calculus which is 
developed in Section 2. Rudiments of such estimates without asymptotic formulae 
but with estimates for remainders appeared in |RS3] . while elements of the calculus 
(without proofs) were announced in [RS5] . In this sense, the present paper is the 
continuation and the development of the subjects of authors' papers |RS3] and [RS4J. 
In particular, the first part of this paper will provide a calculus background to the 
better understanding of estimates of the remainders in [RS3j . this allowing further 
applications to similar problems with certain structures in operators. 

Thus, the first aim of this paper is to develop the calculus necessary for these 
purposes, at the same time investigating the independent question of what are the 
minimal decay assumptions for the global calculus to exist. Consequently, we are 
also interested in global weighted estimates for Fourier integral operators in Sobolev 
spaces, and in applications to smoothing problems. This is especially important in the 
investigation of global smoothing properties for equations with variable coefficients, 
a topic that will be treated in detail elsewhere. 

We will be dealing with operators T that can be globally defined by 

(1.1) Tu{x)= ! I e i ^^-y^a(x,y^)u(y)dy^, 



with a real- valued phase function cf) and a smooth amplitude a(x,y,^), satisfying 
estimates 

\d:d;dja(x,y,0\ < C a ^(x) mi (yj m2 (O m3 , for all x,y,£e E", 

for all multi-indices a, (3, 7. Here and in the sequel we use the notation <?£ = (27r)~™<i£ 
and (x) = (1 + |x| 2 ) 1//2 . As usual, we understand these operators as limits in <S(IR n ) 
as e — > 0, of operators with amplitudes a e (x,y,!;) = a(x, y,£)7( e 2/> e 0> where 7 G 
C^°(M. n x M") equals 1 near the origin. These operators can be naturally extended to 
iS'(IR n ). We will adopt the minus sign in the exponent for some technical convenience, 
and also in order to make the exposition directly applicable to pseudo-differential 
operators, thus providing new weighted Sobolev estimates for pseudo-differential op- 
erators as well. 

Let us note that operators with amplitude a(x,y,^) = a(x, £) independent of y, 
and with representation f 1 1.11) locally in x, appear as solutions to the Cauchy problem 
for strictly hyperbolic equations. Their properties have been studied extensively (see, 
e.g. Kumano-go |Kulj . and also Stein [SU Chapter IX]). A global in x representation 



with amplitude satisfying some decay has been also used to study global properties 
of solutions to the Cauchy problem. Thus, problems with phases and amplitudes 
in so-called SG-classes were investigated by Coriasco [Cor] . In jBBRj. Boggiatto, 
Buzano and Rodino investigated global problems with phases and amplitudes from 
some polynomial classes associated to the Newton polygon of the hyperbolic opera- 
tor. In comparison with these results, we will impose much less restrictive conditions 
on phases and amplitudes. In addition, from our theorems on compositions, it will 
be easy to see that when the original amplitude satisfies some additional decay con- 
ditions, the amplitude of the composition will satisfy the same conditions too. In 
this way we will automatically recover global calculi developed for SG-operators by 
Coriasco [Corj . for operators with certain polynomial weights by Boggiatto, Buzano, 
and Rodino [BBR] . and to a certain extent for operators of Gelfand-Shilov type by 
Cappiello |Caj . as special cases of our calculus. Below we will give an example of 
such an argument. 

A special case of operator T occurs when <p(x, £) = x ■ £. In this case T is a 
pseudo-differential operator. In Corollary 12.101 we will show that if derivatives of the 
amplitude a(x, y, £) decay in the second variable, T can be defined by an amplitude 
p(x, £) of order mi + m 2 in x, and m 3 in £, in the usual form of pseudo-differential 
operators. 

To model weighted Sobolev spaces, we will use a convenient class of operators, the 
so-called SG pseudo-differential operators. We will say that p = p(x, £) G SG tut2 (W l x 
R") if p is smooth and if for all a, (3 we have the estimate 

\d^p(x,0\ < C aP (x) h - m (0 t2 - la \ for all x,£ G W\ 
The corresponding pseudo-differential operator P G Op(SG tl)t2 ) is defined by 



1.2) Pu(x)= [ [ e«*-*>tp{x,ZMy)dy<%. 



Note that in principle one can consider more general SG-amplitudes a(x, y, £) instead 
of p(x, £) in (11.21) . i.e. smooth functions a, for some mi,m2,m 3 G M satisfying 

\d^ x d]a(x,y,0\ < CV 7 <zr H %r H7l <£r 3 - H , for all x,y,Z G R n , 

and all multi- indices a, j3, 7. In |Co2j . Cordes shows that such symbols can be reduced 
to p — p(x,£) G SG Sl+S2,S3 with an asymptotic formula 

Let us show that we automatically recover such (and other similar) results in the 
framework of our calculus. Indeed, from the asymptotic formula of Corollary 12.91 
we see that if the amplitude a(x, y, £) is an SG-amplitude, then p(x, £) must be 
an SG-amplitude as well. This type of argument can be used with different decay 
assumptions (such as polynomial, quasi- homogeneous, Gelfand-Shilov, etc.), to yield 
similar conclusions. 

We note also that the adjoint operator T* of T from (11.11) is given by 



T*v(x)= / e^<-^)a(y,x,Ov(y)dydt, 



and it is operators of this type that appear in global smoothing problems ( |RS2j - 

In Section 2 we will develop elements of the global calculus dealing with global 
compositions of operators ( 11.11) with pseudo-differential operators. In the proofs of 
Theorem 12.11 and Lemma [2~3| we will assume that all integrations by parts are carried 
out even number of times, so that we do not have to keep track of the sign in front of 
the integral. Since we make a choice to represent the phase not symmetrically with 
respect to x and y, we will have non-symmetric statements also in our theorems. 
Thus, we will have two different versions of composition formulae for T o P based 
on whether we prefer to keep an arbitrary phase (as in Theorem 12. ip or to have 
the amplitude dependent on less variables (as in Theorem 12. 8p . For the composition 
P oT pseudo-differential operator P interacts with the phase of T, and we will have 
Theorem 12.51 dealing with this case. In terms of the growth properties of the phase 
function, we will usually think of a typical situation for Fourier integral operators 
when the phase <fi is of order one in both x and £ and the growth may be cancelled 
once we take a derivative (for precise assumptions see the corresponding theorems). 

In Section 3 we will apply the established calculus to derive weighted Sobolev L 2 
estimates for our classes of Fourier integral operators. This will extend weighted L 2 
estimates established by the authors in [RS4j to weighted Sobolev spaces over L 2 , the 
setting important for smoothing problems. 

In Section 4 we will apply the calculus and weighted estimates to establish a number 
of new smoothing estimates for hyperbolic equations, relating time and space decay 
of solution to that of the Cauchy data. Smoothing estimates for dispersive equations 
have been thoroughly analysed in the literature (e.g. [BD] . [BK] . |CS] . |KY] . [RS6J, 
[51] . Sj , [V], [Wa] . and many other papers). The idea is that integration of the 



solution to an evolution equation with respect to time brings an improvement in the 
regularity with respect to x. This is a phenomenon similar to the one exhibited in 
Strichartz estimates where one observes an improved behaviour in the scale of L p 
spaces. Smoothing estimates are centred at the improvement in Sobolev regularity in 
L 2 , and this improvement depends on the order of the evolution equation. Moreover, 
for applications to the time well-posedness problems for nonlinear equations, one is 
looking for global smoothing estimates, in which one inserts additional weights to 
make estimates work globally in space-time. In this respect, the class of hyperbolic 
equations can be viewed as the endpoint case for smoothing estimates, in the sense 
that there is no improvement in regularity, but a question of the weight remains. In 
fact, the problem of finding the weight can be actually viewed as the most important 
here, since the method developed by the authors in [RS6j allows to extend such 
estimates further to other evolution equations (e.g. to Schrodinger, KdV, and others), 
automatically gaining the improvement in the regularity with respect to x. Estimates 
presented in Theorem 14.11 are a new type of estimates for hyperbolic equations that 
we can derive using the global calculus developed in the first part of the paper. 

If we do not specify the domain of integration, it will always mean that the inte- 
gration is performed over the whole space M™. As usual, we will denote D a = i~^d a . 
We also denote <f£ = (2n)~ n d£ and (x) = (1 + Ix) 2 ) 1 / 2 . Constants will be usually 
denoted by letter C, and their values may differ on different occasions. 



2. Composition with pseudo-differential operators 

In this section we will study the composition of T with pseudo-differential opera- 
tors. Let T be given by f 1 1.1 1) and let P be given by (jl.2p . The composition PT is 
then of the form 

(PTu)(x)= f f e^ x - v)v p(x,r])(Tu)(y)dydri 

JR n it" 

(2.1) = f f I f e l ^ x - y) ^ + ' i> ^ ) - z ^ P {x,7 1 )a{y,z,i)u{z)dzdidyd7 1 

JR n JR n JR n JR n 
JR n JR n 

where 

(2.2) Cl {x,z,£)= f [ e iWy ^-^ x ^ +(x ~ y) - v) a(y,z,^)p(x,r])dydrj. 

JR n JR n 

The composition TP is of the form 

{TPu)(x) = f f e^ x ^-yVa(x,y,0(Pu)(y)dydt 

JR n Jl» 

(2.3) = f f f f e lWx ^- y<+{v - z)v) a(x,y,Op(y,ri)u(z)dzdr]dydC 

JR n JR n JR n JR n 

= f f e iWx ^~ z - ,l) c 2 (x,z,ri)u(z)dzdri, 

JR n JR n 

where 

(2.4) c 2 (x,z,ti) = c 2 (x,ti)= I I e i{ ^ x ^-^ x ^ +y ^-^a{x,y,Op{y^)dydi. 



n ./TSn 



n .nan 



Note that this representation of TP makes the amplitude dependent on two variables 
only (independent of z). This will be used later in Corollary 12.91 We also note 
that both amplitudes C\ and c 2 are of the same type, if we interchange variables x, y 
with variables £,77, respectively. The difference is that the integration is carried out 
with respect to one or two variables of the amplitude a in formulae for c\ and c 2 , 
respectively. We will use this fact and will treat both cases simultaneously in Lemma 

IO 

Let us now use another representation of TP to obtain a slightly different estimate 
for its amplitude. In fact, one can readily see that we can represent TP as 



'i J lttn 



(TPu)(x)= I I e i{ ^ x ^- z<) c(x,z,i)u(z)dzdi, 
with amplitude 
(2.5) c(x,z,£)= I I e i to-*M"-Va{x,y,Z)p(y i Ti)dy<lTi. 



n ./IB" 



Note that in comparison with representation (12.3jl - (12.4p . the amplitude here depends 
on three variables, but there is no entry of the phase. This allows to treat more 
general phases and amplitudes, leading to the following: 



Theorem 2.1 (Composition TP). Let operator T be defined by 
Tu(x)= [ [ e i ^ x ^-y^a(x,y^)u(y)dy^. 



n .linn 



Let the phase <p — 0(^0 be any function and a = a(x, ?/,£) E C°°(M. n x R™ x 
satisfy 

(2.6) \d:d;dJa(x,y,0\<C a ^(x) mi (y) m2 (O m3 , for all x,y,£ ER n , 
for all a, /3, 7. Let p = p(y, 77) E C°°(M. n x R") for all a, (3 satisfy the estimate 

(2.7) \d^p(y, V )\ < C^iy^iv)^ 1 , for ally, 7} E R". 
Then the composition B = T o P(x, D) is an operator of the form 

{TPu){x) = [ [ e iWx ^- z<) c(x,z^)u(z)dzGt^ 



n nx 10>ri\ 



n ./IB" 



with amplitude c(x, z, £) satisfying 

(2.8) \d^ z djc(x } z,0\ < C afh (xr(z) m ^(O m3+t \ 

for all a, /3, 7 and all x,z,£ E R n Moreover, we have the asymptotic expansion, im- 
proving in £; 

c(x,z,£) ~2^ — j-9^ [a(x,y,Qd%p(y,£)] \ y=z . 

a 

Definition 2.1 (Improving asymptotic expansions). In this paper, when we say 
that the asymptotic expansion a ~ YllLi a j * s improving in £, it will always mean 
that for every a, /3, 7 and M there exist N and C such that 

8£$8l (a - JTaA < C > (x) m ^y) m * (^T^ \ for all x,y,£e W 1 . 

Similarly, we can define expansions improving in x or in y. 

Remark 2.2. It is clear that if a has additional decay with respect to some variables, 
so does the new amplitude. For example, if we assume stronger estimate 

\d^d^a(x,y,0\ < C aPl (x) m ^ a \y) m *(Z) m \ for all x,y,£ E R" 
instead of (12.61) . then we can improve estimate (12. 8p as 

(2.9) \d^ z djc(x,z,0\ < C af s,(xr-^(z} m ^(0 ms+t2 . 
Similarly, if we assume 

(2.10) \^ y d]a(x,y,0\ < C a p,(x) m ^y) m ^(aT\ for all x,y,£ e R n 
and 

(2.11) \d^p(y, V )\ < C^y) 11 ^^^ for all y, V E R" 
instead of estimates ( 12.61) and ( 12. 7ft . then we can improve estimate ( 12.81) as 

\d^ z d]c(x,z,0\ < C aPy (xr(z) m ^-^(O m3+t2 . 



In fact, from the expression (I2.5p . we obtain 



and the conclusion (12.91) follows from (12.81) if we regard (x) a (d^a)(x, y, £) as a(x, y, £) 
in (12.61) . The same argument is true for the case of assumptions ( 12.101) and (12.111) 
since we obtain 



d g .c(x,z,£)= / / e i to-*>to-Q{d v .a)(x,y,Z)p{y,v)dyitv 
Zjc{x,z,£)= \ I e l{y ' zy{ri ' y j a(x,y,^)p(y,r])dydr] 



" ./IB" 



+ / / e^-^-^a(x,y,0(d rij p)(y,v)dyd V , 



mi 

n ./Wn 



from the expression (12. 5p again. From this we can conclude estimate (I2.9P similar to 
the argument before. 

Proof of Theorem\2Jl Let \ = x{v,z) e C°°(lR n x ]R n ) be such that x(v,z) = l for 
\y — z\ < k/2, and x{Vi z) = for \y — z\ > k, for some small < k < 1 to be chosen 
later. Using this cut-off function, we decompose c = c 1 + c 11 , where c 1 and c 11 are 
of the form (12.51) . but with amplitude a replaced by amplitudes a 1 = (1 — x) a an d 
a 11 = x®, respectively. 

1. Estimate for \y — z\ > fc. We integrate by parts with transposes of operators 
t L T1 = —\y — z\~ 2 A v and t L y = (£)~ (1 — A y ), for which we have 

Thus, we obtain 

c 7 (x,z,0= [[e-tb-KL* [e^- z ^a\x )y) i)L^ V {y )V )\dyd V 

e i(y-z>(v-O e -i(y-zy« L N2 [e^v-^a 1 {x,y,i)L^p{y,r])\ dyd V 
iz . (v _ c) {l-A v ) 3 

< [ e ^-(.-O e -i(,-)-. L iV 2 [e 4 ^V(x,y,0O(^)]]^- 



Now, if we use that after all differentiations, e l ^ y z '"° in the last integral cancels, we 
get 

|c 7 0r,z,OI < c ff (z)- 2Ni (y) 2N3 (0~ 2N2 (v) 2N2 ^) mi (yr 2 (O m3 x 

J J \y—z\>k 

x\y-z\-™i(vt(rjf>- Wx dvdri 

r I \2N 3 +m2+ti r 

< C(x)^(z)- 2N3 (O m3 ' 2N2 I ,., , 2jV dy J (r]) 2N2+t2 - 2N ^ V 



mi i ^\-2N 3 i C \m 3 -2N 2 / \yj_ j„, / / ^\2N 2 +t 2 -2N 1 _ 

L, _ ?\2N 1 
y-z\>k \y /0 \ 

(y + zf 
y \>k \y\ 

.mi / \m 2 +ti / j-\m3-2N 2 / / „,\ 2N :i +m 2 +t 1 -2N 1 



.2Af 3 +m 2 +ti 

< c( X r(z)-^(o m3 - 2N2 1 tfT f„ ffl — ^ 



<C(x) mi (z) m2+tl (£) m3 2 / { y y m+m2+tl -' m dy, 

J\y\>k 

where we have used the estimate (y + z) < C(y)(z) for ./V3 such that 2Ns+rri2+ti > 0. 
This gives the desired estimate for c 1 for JVi ^> Ns,N 2 , and ./V2 > "023. 

For the estimate of derivatives of c ! (x, z,£) we first note that derivatives with 
respect to x do not have any effect on our estimates. Derivatives with respect to 
z introduce the factor of 77 — £, which is taken care of by choosing large Ni,N 2 . 
Derivatives with respect to £ may introduce the factor of y — z, which is taken care 
of by taking large N±. 

2. Estimate for \y — z\ < k. We set 77 = £ + 6, so the Taylor expansion of p with 
respect to 77 gives 

p(y^ + 6)= J2 6 ^p(y,0+ E CJ a r a (y,i,6), 

\a\<M ' \a\=M 

1 

Moa„ 



r a {y,M)= / (l-t) M d?p(y,Z + t9)dt. 







Then 



c n (x, z, = J J e (y ~ z > e X (y, z)a(x, y, S)p(y, £ + 9)dyd9 

= Yl ~ C a(x,Z,^,)+ J] C a Ra(x,Z,g), 



\a\<M \a\=M 

where 

R a (x, z, = ^ ^~'> 9 r X {y, z)a(x, y, S)r a (y, f , 0)dj/d0, 

c Q (x,z,0 = JJ e^~ z > 6 e a [x(y,z)a(x,y,Z)d%p(y,Z)] dydQ 

= Dy [x(v,z)a.(x,y,g)G%p(y,g)] \ y=z . 
From this we readily obtain 

|c a (o;,«,0l<<a;> mi ^> ma+ ' 1 <0 m8+ft, " |a| - 

It is easy to see that derivatives of c a satisfy the same estimates. For the asymptotic 
formula we note that c^x(y, z)\ y=z = for all a ^ 0, so all derivatives must fall 



on the product a Of p. Finally, it is easy to see that derivatives of c a satisfy similar 
estimates as c a . 

3. Estimate for the remainder R a . Let p G C£°(M n ) be such that < p < 1, 
p(x) = 1 for \x\ < e/2 and p{x) = for \x\ > e, for some e > to be chosen later. 
We decompose 

R a (x, z, = Ri(x, z, + Ri'ix, z, 

e<y~ z > e e a p (JL) X (y, z)a{x, y, Or a (y, £, 9)dyd9 



+ f f e iiy ~ z) - e 6 a (l-p( 



777) )x{y,z)a(x,y,C)r a (y,^,9)dyd9. 



3.1. Estimate for \9\ < e(£). For sufficiently small e > 0, for any < t < 1, 
(£ + t9) and (£) are equivalent, which means that \d^r a (y,^,9\ < C(y) x (£} 2 ~' a ' J for 
\9\ < e(£), and all a,/3. We have 



Ri(x,z,0 = JJe liy - z)e D« 

\Ri(x,z,0\<C(xr(O mi+t2 - H 



P\T7\) X(y,z)a(x,y,£)r a (y,£, 



dyctO, 



d6 / (y) 

mi /j-.\ m 3 +t 2 -|a|+n / \ m 2 +*i 



m 2 +ti 



dy 



<c(x) mi (O m3+2 ~ m+n (zy 

which yields estimate (12.141) for M = \a\ > n and large z. 

3.2. Estimate for \8\ > e(£). Here we will use that fact that for M = \a\ > £2, the 
function r a (y,^,9) and its derivatives can be estimated by (y) 1 . We will integrate 



by parts with operator L y 
have 



\9\- 2 A y , for which we have L N e i{ - y - z ^ B = e^" z K We 



R*(x,z,0 



e ^y 



P f "Try ) ) X(y,z)a(x,y,^)r a (y,^,9) 



dyd9 



,i{y-z)-6 



\o\ 



-2N A N p,a 

y y 



P f -Try ) ) x(l/> z)a(s, 2/, O r «(2/ 5 ^ #) 



dyd9, 



which implies 



i^(x,^,oi < c(x) m ^o mi / ir 2A ^ / (y) m2+ti ^ 

J\e\>e(0 J\y-z\<k 

<C(x) mi (O m3+n ~ 2N (z) m2+tl - 

Derivatives of i?„, R X J with respect to z and £ may introduce additional powers of 9, 
which amounts to taking more ^/-derivatives under the integral and does not change 
the estimates. This completes the proof of Theorem 12.11 □ 

Now we will prove a lemma which will be used in deriving estimates for amplitudes 
in compositions PT and TP in fl2.1l) - fl2.4p . For this purpose, we introduce amplitudes 
that may depend on four sets of variables, and we will be able to choose different sets 
when dealing with different composition formulae in the sequel. 



Lemma 2.3. Let ip = tp(x,£) G C°°(R n x R n ) be such that 

(2.12) Ci(0 < (V^(x,0) < C 2 (0, /or a 7 'x,£ G M n , 
/or some C\,C% > 0, and such that for all \a\, |/3| > 1 we aat>e 

(2.13) |W(*,£)I < Ca<0, |^V(x,e)l < Ca/3, M a// x, e e » n . 
Le£ 6 = b{y, £, z, n) G C°°(R 4n ) satisfy 

\8$dld!a$b{v,Z,z } TJ)\ < c aM s{y) mi {i) m2 {z) m ^ v r^ s \ 

for all a, /3, 7, S, and all y, £, z, 77 G M n . Le£ g = g(x, 77) G C°°{R n x E n ) /or a// a, (3 
satisfy 

\d«dPq(x,T})\ < C a p(x) h (ri) t2 ' m , for all x,r] G R n . 
Then the function 

c(x 1 z 1 £)= f f e l ^ y ^-^ x ^ + ^- y ^b(y,^z,ri)q(x,ri)dy(!7 1 

satisfies 

(2.i4) \d:d^ z djc(x,z,o\ < c a ^(xr +h (zr*(o m2+t2+m \ 

for all a, /3,7, and a// x, z, £ G R n . Moreover, we have the asymptotic expansion, 

improving in £: 

(2.15) 

,— (|a|+|/9|) 

c{z,z,£) ~ £ ^-^ a g(x, V^(^0)^ Q+/3 [e^*^%,£,z, V^(*,0)] U, 

wnere *(&, y, = V(l/, - ^(s, + (a; - y) ■ V^ar, f )• 

Proo/. Let x = x(&,2/) e C°°(R n x R n ) be such that < x < 1, xfay) = 1 for 
\x — y\ < k/2 and xi x i u) — f° r 1^ — 2/1 > &• Let 

c(z, «, = c 7 (x, z, f ) + c /7 (x, z, 

+ J J e my£)~^)H*-y)-v) x ^ y)b{yj e> Zj r/)g(X) ^ % 

1. Estimate for \x — y\ > fe. Here we will integrate by parts with two operators, 
t t ■ \ v "^o — Vi ' o t t — v 

v = ~ l U^^ 2 v= (VyiKv,Z)) 2 -iW(v,ZY 

so that we have L^e 1 ^'^ 1 = e^"^, L^e ilp{y & = e nKv.0. Note that in view of our 
assumption (12.121) on ip, we have 

\(V y iP(y^)) 2 -iA y i;(y,0\ > |(V^Q/,0>| 2 > C^f ■ 
Therefore, 

c?(x,z,£) = jjS+M)-4>M)+*v)) L N* (e-^L^ [(1 - X (x,n))6(n,e,z,n)g(x,n)]) dydrj. 
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Further, we integrate by parts with operator (1 + \x\ 2 ) 1 (1 — A,,), for which (1 + 
|x| 2 ) _1 (l - A^e"" 7 = e ix -v. This gives 

<?{x, z, = If e ™-v e iMy,Z)-^)) ( x )- 2N *(i - A V ) N3 L^ 
(e~*™L* [(1 - X (x, y))b(y, £, z, V )q(x, V )]) dyd V . 

Since Ni ^-derivatives of b(y,£,z,r))q(x,r)) in view of our assumptions give an im- 
provement in i] by (r/)~ 1 , we get 
(2.16) 

1^,01 < // (x)- 2N3 (y) 2N3 (0- 2N2 (v) 2N2 

J J\x-y\>k 

\x - y\' 2Nl (y) mi (0 m2 (z) m:i (*>* (v) m4+t2 ' Nl dyd v 

< C {z) m *(x)^- 2N *{t;T*- 2N * / ( V ) 2N2+m * +t2 - m d V / \ y) _ dy 

J J\x-y\>k \ X y\ 

< c(z} m3 ( X ) h - 2N3 (o m2 ~ 2N2 / {y | L — */ 

<C(z) m3 (x) tl+mi (O m2 ~ 2N2 i 

where we used (y + x) < C(y)(x) with N 3 such that mi + 2N 3 > 0, and 2Ni > 
rrii + 2-/V 3 + n, 2N 2 + m 4 + t 2 + n, to make both integrals finite. 

Derivatives with respect to z do not change estimates. Derivatives with respect to 
x and £ may introduce factors estimated by powers of (£), (77) and \y — x\, which can 
be taken care of by choosing large iVi, JV2. 

2. Estimate for |x — y\ < k. We will make a Taylor expansion of q{x,n) and 
%, £, 2, V) in ^ = V x ^(a;, + 0: 

<?M = E ^g(x,V^(x,0)+ E C<#V a (z,£,0), 

a|<Mi ' \a\=Mi 

(2.17) l/3|<^2 ^' |/3|=M 2 

r a (x,£,9)= I (l-t) M id«q(x,VMx,0+t6)dt, 
Jo 

s p (y,Z,z,9)= I (l-t) M2 dlb(y^,z,V x ^,0+te)dt. 
Jo 



Let 



Then 



*(x, y, = V(y, - V>(z, + (x - y) ■ V x ^(x, £)■ 
c /7 (a:,^,0= E "Tgi ^(^ z ' + E CaC^-R^^,^,^)^ 

|a|<Mi,|/?|<M 2 '^' |a|=Mi,|/3|=M 2 
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where 
Rceix, z, = J! e^-y> e e m ^e a+ e X (x, y)sp(y, £, «, 0)r a (x, £, *)«*!/*», 

c af3 (x,z,£) = He^-y> e 9 a+ ?x 

[e^ x ^ X (x, y)d%b(y, £, 2, V^fo £))^ g (*, V^(rr, 0)] ^ 



x l f ' • ^ - 

1 2/=a; • 

Since y = x is the zero of order two of ty(x,y,£), it follows that V ye l ^^ x ' y ^'\ y=x = 0. 

M + I0I 



Inductively, it follows from (12. 13[) that |92e l **' a: ' J/ '^L =a ;| < C(£)' 7 '' . Hence we obtain 



\Cafi(x, Z, | < C(X) 41 (V,V>(*, fl)*"'"' (0^~ (^) mi (0" 12 fcP <V.^(X, O)" 14 "" 31 



< a(a;} mi+ * 1 (^) m3 (0 m2+m4+ * 2_ ^ M . 



where we have used that (V x ip(x,^)) ~ (£) by (I2.12p . Derivatives of c a p with respect 
to x or £ do not change estimates when applied to q or b by assumptions of the lemma. 
When applied to dy + Pe iq, ( x,y '^\ y=x they do not change estimates since \ded^ip(x, £)| < 

C a j3- 

3. Estimates for the remainder R a p. Let p G C£°(K n ) be such that p(x) = 1 for 
\x\ < e/2 and p(x) = for |x| > e, for some small e > to be chosen later. We 
decompose 

(2.18) 

Rafiix, Z, £) = i?^(x, Z, + i£ J (x, Z, f ) 

= || e M '<V (|y) D^ [ X (x, y)s p (y, £, *, 9)r a (x, £, (9)] dydd 

+ 1| e «*-v)* (\ _ p (o\ \ uc+p [x{X} y)s ^y } £ } Zi e)ra{x7 ^ 6 )} dyde. 

3.1. Estimate for \6\ < e(£). For sufficiently small e > 0, for any < t < 1, 
(V x , ?/'(x,£) + t6) and (£) are equivalent. Indeed, if we use the inequality 



(x) <l + \x\ < V2(x), 
we get 

(VMx,0+t9)<(C 2 V2 + e)(0 

V2(vMx,0 + to) >i + |v^| - \e\ > (v.v) - |0| > (Ci - e)(0, 

so we will take e < C\. This equivalence means that for \9\ < e(£), functions 
r a (x,^,9), S/3 (?/,£, 2,0) and all of their derivatives are dominated by (x) x (£) 2 ~' Q| 
and (i/)" 11 ^ m2+m4 ~iPi ( z ) m3 ^ respectively. We will need two auxiliary estimates. The 
first estimate 

12 



d] 1 (p (|y) r a (x,e,fl)) < CJ2 %P (Jfi) ^r^x^O) 



(2.19) 



<5<7 

5<7i 

< C(a;> tl <£> ta ~ |o ' H711 



l<5|/M*2-|a]-|7l-<5| 



follows from the properties of r a . Before we state the second estimate, let us analyse 
the structure of d^e l ^^ x,v ^' . It has at most |a| powers of terms V y ip(y, £) — V x ip(x, £), 
possibly also multiplied by at most |a| higher order derivatives dyip(y, £), which can be 
estimated by (\y — e|(0) ^ (12. 131) . The term containing no difference V y ip(y,£) — 
V x i/j(x,£) is the product of at most |a|/2 terms of the type d^iftty^), which can be 

estimated by (£) in view of ( I2.13p . Altogether, we obtain the estimate 

|^ e i *^«|<c a (i + |j/-zK0) |a| <0 w/2 - 

The second auxiliary estimate now is 

(2.20) <(i + \y-x\(0) laW K0^(y) mi (0 m2 (z) m H0 m4 ~ lPHj21 

<(i + |y-x|(0) H+l/3| (y) mi ^) m3 (O m2+m4+M ^ M_l721 - 

Now we are ready to prove a necessary estimate for R^Jx, z, £). Let 



L« 



(i - (0%) 



j^N e i(x-y)-e _ e i(x-y)- 



Integrations by parts with L e yield 



#iaO>^0 



D i(x-y)-9 



;i - (o% 



\A/ 



i + <0>-y| a r 

E„ /t\l7l|+l72| 

°7i72\S/ 



J(x-y)-6 



(i + (0 k - i/l a ) 



some *yi ,72 : 
| 7 i| + | 7 2t<2Af 



D 



a+/3 



X(x, y)d 1 e 2 sp(y 1 £, z, 6)df [ p ij-r ) r a (x,£, 9) 



( 



> Gfa/c?#. 
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Using estimates (j2.19p . (12.201) for r a , Sp, and the fact that the size of the support of 
p(9/(0) i n ^ i s estimated by (e(^)) 71 , we obtain the estimate 



l+|7i[+|7a|/„\tl /c\*2-|Q|~|7i|/^\ m 3/^\"t2+m 4 + |Q| l/3| -|7 2 | 



hi\ + \l2\<2N 

'i«-i<* (i + (0 2 \x - y\ 2 ) N 



< C 7 r\' 1+mi /c\ t2 + m2 + m - 4 + 2re ~ 2 /?A m3 



if we choose 2iV > M1 + M2+TI. In the last estimate we have used that for \x — y\ < k, 
(x) and (y) are equivalent for large x. 

Derivatives of R^Jx, z, £) with respect to x and £ give an extra power of 9 under 
the integral. This amounts to taking more ^/-derivatives, giving a higher power of (£). 
However, this is not a problem if for the estimate for a given number of derivatives 
of the remainder R^Jx, z, £), we choose M x = |a| and M 2 = |/3| sufficiently large. 

3.2. Estimate for \9\ > e(£). Let us define 

w(x, y, £, 0) = (x - y) ■ + *(x, y, = (x - y) ■ (V^(x, £) + 9) + i/;(y, - if>(x, £)■ 

From (12.121) and (12.131) we obviously have 

| V y oo\ = \-9 + V y ij- V x ^| < 2C 2 (\9\ + (£)), 
(2 ' 21) IV^I > \9\ - \W y ,p - V^l > ^1 + (| " Col* - y|) (0 > C(\9\ + (£», 

if we choose fc < e/2C , since |x — y| < fc in the support of x- Let us denote 
a 7l (x,y,£) := e -<*(*.v.O£)TV*(*.v.O x 72 (*,y) := ^ y 72 x(*,y), ^ 73 := DJV 

For any z/ we have an estimate |<9y<r 7l (x,y,£)| < (£) , because of our assumption 
fl2~T3l that \d%ij)(y,g)\ < C„(f). For M x = |a| > t 2 and M 2 = \p\ > m 4 , we also 
observe that 

(2.22) \r a {x,{,e)\ < C a (x) l \ |^ 73 (y,£, z,^)| < C $ {yf* (&»{*)"*' 

Let us take t L y = i\ V y u}\~ 2 Y^j=i(^yj u )^yj- ^ can ^ e shown by induction that operator 
Ly has the form 



(2.23) L N y = — 1— £ P^, P„ iJV = J] c^ (V^)^ w • • • fl{"u; 



a ~ ' |i^|<iV |^|=2iV 
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where |/x| = 2N, \5A > l,J2i \ 5 j\ + W\ = 2N - lt follows from (12331) that \P U>N \ < 
C{\9\ + (C)) ZN ■ By Leibnitz formula we have 



R I J (x,z,O = 

e i(x ~ y) - 6 I I .- ,>[ — ) )r n U.Z,.0)x 



*/9 



(2.24) 



x d^ [j*^ x (x,y)sp(y,t;,z,e)] d y m 

7i+72+73=a+/9 

x 5Z d y( (T 'ji( x ,y,Oxi 2 (x,y)sp'y 3 (y,£,z,&))dydO. 

7i+72+73=a+/3 

It follows now from (12.221) that 

I35(*,*,fll < c I f {\e\ + (0)- Ar (^) il (0 |Q|+l/3| (i/) mi (O m2 (^r 3 ^^ 



•/|0|>e(£W|j/-x|<& 

<c(x)' 1+mi (^) m3 (O m2+|a|+l/31 / \e\- N de 

< C(x) tl+mi (z)" 13 (c\ m2+ \ a \ + \< a \ +n ~ N 

which yields the desired estimate for iV > m 2 + \a\ + |/3| + n. For the derivatives of 
R T JJx, 2, £), similar to Part 3.1 for R^n, we can get an extra power of 9, which can 
be taken care of by choosing large N and using the fact that \x — y\ < k. The proof 
of Lemma 12.31 is now complete. □ 

In the proof of the second part of Theorem 13.51 we will need to have stronger decay 
properties for the amplitudes, so now we will formulate the corresponding counterpart 
of Lemma 12.31 

Proposition 2.4. Under conditions of Lemma \2.3\ assume that stronger estimates 

(2.25) \d^d:d s ri b(y,^z, V )\ < C a ^s(y) mi - lal (0 m2 (zr 3 (v) m4 - lSl 
and 

(2.26) \9*tfq(x,Ti)\ < C aP (xf- lal ( V ^- m 

hold. Assume also that phase function ip(x,C,) satisfies in addition estimates 

(2.27) |W(x,0l < Caix) 1 -™®, \d:d^(x,0\ < C aP (x)^ al , 

for all multi-indices a and (3 with \/3\ > 1, and all x, £ G MJ 1 . Then we can improve 
estimate (12.141) as 

(2.28) \^djc(x,z^)\ < C a/3l (x) mi+tl - lal (z) ms (0 m2+t2+m4 - 
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Proof of Proposition pT^j In the case of proving estimate (12.281) . we have to go back 



and modify the proof of Lemma 12.31 to adopt it to the required decay properties. We 
will assume the notation of the proof of Lemma 12.31 and will indicate the differences 
for the present case. The main difference is that now we will work with a different 
cut-off function x = x{ x ->v) e C°°(IR n x W 1 ) such that < x < 1, x{ x tV) — 1 f° r 
| a; — y\ < |(x) and x( x iV) = f° r \ x ~ y\ > k(x), for some small < k < 1 to be 
chosen later. 

1. Estimate for \x — y\ > k(x). The proof goes similar to the corresponding 
proof of the first part of Lemma 12.31 In the estimate for c 1 in (I2.16P we observe 
that the domain of the integration in the last integral changes from \y\ > k to 
\y\ > k(x), thus introducing the additional factor of (x)~ for any N, if we choose 
2Ni > mi + 2N3 + n + N. Therefore, the part c 1 of the amplitude c is actually 
decaying as any negative power of (x) in this case. 

2. Estimates for [re — 3/j < k(x). In Part 2 of the proof of Lemma [2.3l we established 
the formula for the terms in the asymptotic expansion of c 11 , given by 



z a p(x,z,0 = d£q(x, V^(x,0)^ [e m ^x(x,y)dPb(y,Z,z, V*V(*,0)] 



y=x- 



Now, differentiation of d?q(x, ^7 x ip(x, £)) with respect to x brings the decay by (x)~' 
in view of assumptions (I2.26[) and (I2.27p . Differentiation of x an d b with respect to 
both x and y also gives a factor of (x)~ on the support of x- Finally, since 

V(x, y } = <P(y, - 00, + (x - y) ■ V x <j>(x, 0, 

we have V^e^^'^l^^ = and D%V(x,y,£)\ y=x = D2<f>(x,£) for \a\ > 2. Thus, by 
the product rule applied to the derivatives Dy +l3 , function c a p contains derivatives of 
X and b, multiplied by terms of the form D^<p(x, £) coming from the differentiation 
of the exponential e 1 *^*^ at y — x. Differentiation of such terms with respect to 
x also brings the decay of order (x)~ , thus establishing the required decay property 
(12.281) for functions c a p. 

3. Estimates for the remainder R a p. Here again we follow the proof of Part 3 of 
Lemma 12.31 adopting all its notation. 

3.1. In the region \9\ < e(£), the definition of r a in (I2.17P and assumptions (I2.26P 
and (I2.27P imply that d1r a brings an additional decay by (x) compared to that of 
r a . Now, differentiation of R ! a/3 in (I2.18P gives this decay when the ^-derivatives fall 
on r a . If the x-derivatives fall on the exponential e^ x ~ y ^' e , they give additional powers 
of 8, which transform into ^/-derivatives once we integrate by parts with respect to y. 
Since y-derivatives of both x an d sp bring decay in y, and hence also in (x) in view 
of (x) < C(y) on the support of x ( we choose k to be sufficiently small), we establish 
(12^81) for functions R^. 

3.2. In the region \9\ > e(£), we note first that the estimate \V v uj\ > C(\9\ + (0) 
in (I2.2ip is still valid in view of (I2.27p . if we choose k > to be sufficiently small 
compared to e. 

Now, we will show that under assumptions of Proposition \2A\ R^p is actually of 

order (x)~' in x, for any N. We note that we have \dyUj\ < C$(y) (\8\ + (^)), 
so that function P V)N in (I2T23D can be estimated by \P V)N \ < C(y) WhN (\6\ + (£)) 3N ', 
since Yli (1 — l^j'l) = N — (2N — \i/\) = \i/\ — N. On the other hand, functions 
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<7 7l ,X72 an d S/3^ give the improvement by (y) when differentiated v times with 
respect to y. It follows that <9^(<7 7l x 72 S/3 73 ) has growth in y of at most (y) m . 
Consequently, we obtain that the function under the last integral in (I2.24p is of order 
{y) ' +( - m ' A ~^> = (y) m3 ~ J in y. Since y is equivalent to x on the support of x f° r 

small k, we see from ( 12.241) that R^* satisfies the estimate 

(2.29) |i22,(s,z,OI <C(x) tl+mi - N (z) m3 (O m2+lal+m+n ~ N - 

Now, let us show that differentiations (especially with respect to x) guarantee orders 
of other variables as in (12.281) . The definition of r a in (12.171) and assumptions (I2.26J) 
and (12.271) imply that d].r a introduces an additional factor bounded by (x) (£) 
compared to that of r a . Differentiating RfL in (I2.24p with respect to x, by the product 
rule we get that if the derivatives fall on either r a or x, we get the additional decay in 
(x) times possible growth in (9). If they fall on the exponential e^ x ~ y ^ 9 , we integrate 
by parts with respect to y to get the decay in (y) times possible (polynomial) growth 
in (9). Since (y) > C(x) on the support of x, in both cases we get the improvement in 
(x) compensated by possible growth of the same order in (9). Derivatives of e i *( x ' v '® 
may give growth in x, y or £, but this is compensated by the arbitrarily fast decay in 
these variables. Altogether, since in (I2.29J) and in (12.241) we get the overall decay of 
R t Jq to be of any desired order in all x, £ and 9, this compensates possible additional 
growth that we obtain by differentiating R^g. This completes the proof of (I2.28p . □ 

Applying Lemma [2.31 to (12.21) and (12. 4p . we obtain two composition theorems. 

Theorem 2.5 (Composition PT). Let operator T be defined by 

Tu(x)= [ [ e^ x ^-y^a(x,y,Ou(y)dydt 



Let the phase = 0(x, £) G C a 

Ci(0 < (V..0(x,O) < C 2 (0, for allx^ G W 1 , 
for some C\,Ci > 0, and such that for all \a\, |/3| > 1 we have 

\d:cf)(x,0\ < C a (0, \d^<P(x,0\ < C a p, for allx,Z G R n . 
Let a = a(x,y,£) G C°°(M n x W 1 x R") satisfy 

(2.30) \d:d^a(x,y,0\ < C a ^{x) m '{ y ) m %i) m \ 

for all a, (5, 7, and all x, y, £ G M". Let p = p(x, £) G C^R™ x R") for all a, (3 satisfy 

(2.31) \8%d%p(x,£)\ < C Q/3 (x) 4l (O t2_l/31 , for allx,£ G R n . 
Then the composition B = P(x, D) oT is an operator of the form 

Bu{x)= I e i{4,{x ^- z<) c(x,z,^)u(z)dzd^ 



with amplitude c(x, z, £) satisfying 

(2.32) \d:d^ z djc(x,z,0\ < C a ^(xr +h (z) m *(O m3+t \ 
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for all a, (3, 7, and all x,z,£ G M. n . Moreover, we have the asymptotic expansion, 
improving in £; 

w/iere #(&, y, = 0(y, - 0(x, £) + (x - y) • V x (f){x, f )• 

Proof of Theorem \2.5\ Formula for c(x,z,£) is given in (12.11) and (12. 2p . The rest 
follows from Lemma 12.31 with if> = <f>, q = p, and b(y , £ , z , rj) = a(y,z,£,), with a 
different notation for the orders. Asymptotic expansion follows from (I2.15P if we 
notice that we only have terms with (3 = 0. □ 

Remark 2.6. As in Remark 12.21 an additional decay of a implies the same decay 
property of new amplitude. For example, if we assume stronger estimate 

\d%d^a(x,y,0\ < C Q ^(xr(y} m2 -W(O m3 
instead of (I2.30p . then we can improve estimate (12.321) as 

(2.33) \d^dj C (x,z,o\ < c a ^(xr +h (zr- m (o m3+t2 - 

In fact, from the expression (J272]), we obtain 

(zfd^c{x,z^) = f f e iWy ^-^ x ^ +{x - y) ^(z) p (d^a)(y,z^)p(x,ri)dyatri 



■n ./IB" 



and the conclusion (I2.33J) follows from (I2.32J) if we regard (y) (d!jja)(x, y, £) as a(x, y, £) 
in (123UD . 

For the additional decay properties with respect to x we have the following: 

Remark 2.7. Under conditions of Theorem 12.51 assume that stronger estimates 

\aSdSdIa{z,v,t)\ < C aPl {xr^\y) m ^r 



and 



\d:d^p(x,o\<c aP (x) h - H (o t2 - m 



hold instead of (12.301) and (12.31[) . Assume also that phase function <f>(x,£) satisfies 
in addition estimates 



IW(*,0| < c a {x) x - m o, \d:d^(x,o\ < c a0 (xy 

for all multi-indices a and (3 with \/3\ > 1, and all x, £ G M n . Then we can improve 
estimate ( 12. 32ft as 

(2.34) \%d!dlc{z,z,Z)\ < C af s,(xr +t ^ al (z} m2 (0 m3+t2 - 

To prove the statement of Remark 12.71 we use Proposition 12.41 with if) = 0, q = p, 
and b(y, £, z, rf) = a(y, z, £). Assumptions of Remark [2TT1 then guarantee that assump- 
tions of Proposition 12.41 are satisfied, so f!2.28|) implies (12.341) . 

The composition formula for TP is given in Theorem 12.11 under less restrictive 
conditions on the phase. However, the following theorem shows that if the amplitude 
a(x,y,£) has decay properties in y, the amplitude of the composition TP can be 
made dependent on two variables only. 
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Theorem 2.8 (Composition TP). Let operator T be defined by 
Tu(x) = [ [ e^^-y^a(x,y,Ou(y)dy^. 



n ./lD>n 



Let the phase 4> = <j>(x,£) E C°°(R n x R") be such that 

(2.35) d(x) < (V^(x,0) < C 2 (x), for all x,£ E R n , 
/or some Ci,C 2 > 0, and such that for all \a\,\/3\ > 1 we /iai>e 

(2.36) |9f0(x,OI < Cf,(x), \d a x d^{x,i)\ < C a0 , for all x,^ E W 1 . 
Let a = a(x,y,£) E C°°(W l x R n x R") satisfy 

(2.37) \^dja(x,y,^\ < C a/h {x) mi {y) m *- m {O™ 3 , 

/or a// a, /3, 7, and all x,y,£ E R n . Let p = p{x, f ) E C°°{R n x R n ) for all a, (3 satisfy 

\d^ V {x,0\<C^(x) tx - H (^, forallx^EW 1 . 
Then the composition B = T o P(x, D) is an operator of the form 

Bu{x)= I e'^tys, £)«(£)#£> 

with amplitude c(x, £) satisfying 

i^ C (x,oi < c Q p(xr +m ^(o m3+t \ 

for all a, /3, 7, and all x,z,£ E R n . Moreover, we have the asymptotic expansion, 
improving in x: 

j-(\a\+\0\) 

c ^) ~ E mr 9 "KV^(^6,6^ [e^'^a(x, V^(*, £),??)] U, 

where Wfa £, x) = <fi(x, f ) - <j>(x, rj) + (rj - f) • V € 0(x, £)■ 
Proof. Let us rename the variables 

x =ri,y =£,£ = a;, 77 =3/. 
Then formula (12.41) becomes 

c 2 (£', x ') = /7 e W £')-W £')+(*' -v'W^y'^ ^ 7? /) g ( x / ; 7/)^'^', 



if we choose ip(x',£') = <fi(£', x'), q(x r , rf) = p(i]',x'), and b(y',£ r ,z,r]') = a(£',r)',y'). 
Then we can apply Lemma [2.31 to complete the proof. □ 

If we apply Theorem 12.81 with p = 1, we can simplify the amplitude of the operator 
T. Operators of such type is a simplified version of Fourier integral operators and the 
local versions are described by, for example, Kumano-Go [Kulj . Stein [St], and many 
other authors. This allows one to make representations for operators as oscillatory 
integrals to depend on fewer variables. 
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Corollary 2.9. Let T be the operator from Theorem \2.S\ with phase <fi satisfying 
Ii2.35\) . IJ2.36]) . and amplitude satisfying {2.31^ . Then T can be written in the form 






with amplitude c(x, £) satisfying 



for all a, (3, and all x,£ EW 1 . Moreover, 

0\ 



l^c(x,0l<C^(x) mi+m2 (0 m3 , 



c (^) ~ E %1 W< [e^^d^y^)]^^] \ v= „ 



where \1/ is as in Theorem \2.8l 

We can apply Corollary l2.9l to pseudo-differential operators to obtain a the standard 
quantisation for pseudo-differential operators with amplitudes having decay in only 
one of the variables. Different quantisations of operators are possible, including 
symbols of the type c(tx + (1 — t)y, £) for < t < 1, with Weyl quantization given by 
t = 1/2. However, for the purposes of this paper and for simplicity we will restrict 
our attention to the case t = 1. Other cases may be done with small modifications. 

Corollary 2.10. Let T be an operator of the form 

Tu(x) = [ [ e« x - v Ka(x,y,Z)u(y)dy<%, 

JM." JR n 

with amplitude a = a(x,y,£) G C 00 ^™ x W 1 x R") satisfying 

\d^ y dja(x, y ,o\ < c aPl {xr{ y r-^^r\ 

for all a, (3, 7, and all x,y,£ G M. n . Then T can be written in the form of a pseudo- 
differential operator 

Tu(x)= [ e fa «c(x,0w(0d£, 

JR n 

with amplitude c(x, £) satisfying 

\d^ C (x } o\<c a ^r i+m2 (o m \ 

for all a, (3, and all x, £ G M n . Moreover, we have the asymptotic expansion 

j-|/3| 
c{x,0 ~^2—^-d^a{x,y,C)\ y=x . 

Note that from all asymptotic expansions it is clear that if a has additional decay 
with respect to some variables, so does the new amplitude (see Remarks I2.2[ 12.61 and 
12.71 for such conclusions). 
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3. Weighted estimates in Sobolev spaces 
Again, we will be studying operators 



(3.1) 



Tu{x) 



SM^-vQafay&uMdyat 



with real-valued phase 4> = <f>(x,£) G C 00 (R ?1 x R n ) and amplitude a = a(x,y,£) G 
C°°(R n x R n x R n ). In |RS4] . we proved weighted L 2 estimates for such operators, 
also giving explicit estimates on the required number of derivatives of phase and am- 
plitude. In Theorem 13.51 we will give weighted Sobolev estimates for these operators. 
We will say that / G H'^(R n ) if / e S'(R n ) and II S1>S J G L 2 (R n ), where U SuS2 
is a pseudo-differential operator with symbol 7r SliS2 (:r, £) = (x) Sl (£) S2 . We can also 
define operators n sijS2 with symbols 7f sl , S2 (y,0 = (y) Sl (0 S2 - Then it can be shown 
(e.g. [Co2] ) that II_ Sl _*„ = II" 1 ... Also, Sobolev space defined bv IL, s „, i.e. the 



-81,-82 



space of all / G S\R n ) with U S1S J G L 2 (R"), coincides with H S1 ' S2 ( 



The strategy for obtaining estimates in weighted Sobolev spaces will be to use 
the calculus developed in the previous section to reduce the problem to global L 2 
estimates for model cases. We will first formulate global L 2 estimates for operators of 
the form (13. ip . In fact, under certain conditions, we will reduce them to operators of 
a simplified form, for which we have the following result, which appeared as a special 
case of Theorem 2.5 in |RS4] : 

Theorem 3.1. Let operator T be defined by 



Tu{x) 



z*M*<>'*»a(y,Z)u(y)dydt. 



<p(y,g) G C°°(R n x R n ) satisfy 
\dstd v dt<p(y,t)\>C>0, 



<C, 



Let the real-valued phase ip 

(3.2) 

and 

|a£%>(i/,0|<Ca 

for all 1 < \a\,\/3\ <2n + 2 and all y,£ G W 1 . Leta = a(y,£) G C°°(R r 

d^a(y,o\<C a0 , 
| < 2n + 1, and all y,£ G R n . Then the operator T is L 2 



d y d$ip(y,g) 



p 



for all \a\, \(3 
and satisfies 



W 1 ) satisfy 



)-bounded, 



\T\ 



L'^L 2 — C SU P 

|a|,|/3|<2n+l 



a?#a(i/,fl 



L°°( 



l n ) 



Remark 3.2. In Theorem 2.5 in [RS4j . we treated the case that conditions on the 
phase (p are imposed only on a subset of R™ x R™, which contains the support of a. 
In such case, however, condition (13.21) must be replaced by 

(3.3) |%>(x,£) -d£<p(y,g)\ > C\x-y\, \d v (p(y,£) - d y (p(y,r])\ > C\£ -rj\. 

We remark that if condition (13. 2p is satisfied for all y, r\ G R™ as in the case of Theorem 
13.11 then the global inverses of V^(',0 and V y (p(y, •) exist, and condition (13. 3p is 
implied by the mean value theorem applied to them. 
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Now, we can use this result to treat operators of the form (13.11) : 
Theorem 3.3. Let T be defined by 

Tu(x)= [ [ e i ^ x ^-y^a{x,v,Ou{y)dydt 



n JVn 



Let the real-valued phase = 0(x, £) G C°°(IR n x IR n ) for some positive constants and 
all \a\, \/3\ > 1 satisfy 

\detd x d^(x,0\ >C >0, 

(3.4) 3x :\dl<l>{xp^)\<C t 



Pi 



\d:d^(x } o\<c a ^ 



for all i,(6l n . Let amplitude a = a(x, y, f ) G C°°(IR n x R n x W 1 ) for some m G R 
satisfy 

\d"d^dja(x,y,^)\ < C a p 1 (x) m (yy 111 ^ , for all a, (3, j, and all x,y,£ el". 

T/ien T zs bounded from L 2 (R n ) to L 2 (R n ). 

Remark 3.4. There is an alternative formulation of Theorem 13.31 if we impose con- 
ditions ( 13.41) on the support of a only, similar to Theorem 2.5 in [RS4] (see Remark 
13.21) . However, in this case we also need to assume that 

(*><C<V^(a;,0> 

on the support of a as well, since it does not automatically follow from conditions 

432D • 

Proof of Theorem \3.3[ By Corollary 12.91 we can reduce T to 



Su(x) = / e^^c(x,0£(0^ 

with amplitude satisfying |<9"<9fc(x,£)| < C Q/ g, for all a, (3. We only need to check 
that conditions (I2.35P of Theorem 12.81 on the phase are fulfilled. For xp as in (13.41) . 
and some z we have 

d^(x,0 - d/*4>{xp,£) = V x d^(z,0(x - xp). 

It follows from ([31} that 

(d^(x,0) < l + |^0(x,OI <l + C\x-xp\+C <C(x). 

So assumption (12.351) is fulfilled if we show that (x) < Ci(Vf0(a;, £)}. From the first 
inequality in (13.41) . we obtain inequality 

|V^(x,0 - V^^/3,01 > C\x - Xp\, for all x,£ G R n 

(see Remark 13.21) . Now estimate (x) < Ci(V^0(x,^)) follows from the second in- 
equality in (13. 4p . The adjoint S* of S satisfies conditions of Theorem 13.11 so it is 
bounded in L 2 (R n ). □ 

Now we will apply this as well as the calculus developed in the previous section to 
obtain estimates in weighted Sobolev spaces. 
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Theorem 3.5. Let operator T be defined by 

Tu(x)= [ [ e iWx ^-y^a(x,y,£)u(y)dyat£. 



n J»n 



Let the real valued phase = 4>(x,£) £ C°°(R n x W n ) for all \a\, |/3| > 1 and all 
x, C, G M" satisfy 

(3.5) |detcW(x,OI > Co > 0, |^(a;,OI < C«<£>, |^0(a;,OI < <V 
Assume one of the following: 

(1) For all a, (3, and 7, assume 

\d:d^dja(x, y ,o\ < c a ^(xr( y r- m (o mi , 

and for all \/3\ > I, assume 

3xp: 1^0^,01 <C P 

for allx,y,£ G M n . 

(2) For all a, f3, and 7, assume 

\d:d^a(x,y,0\ < C a ^(x) m ^ al (yr(O m \ 
and for all a and \/3\ > 1, assume 

\d^(x,0\ < Caix) 1 -™®, \d:d^(x,0\ < C a0 {x) l -W, 
for allx,y } £ G M n . 
Then T is bounded from H ai > a2 (R n ) to #*i-«u-«»a.-a-«»s (R n ) , for all s u s 2 G R n . 

Proof. Let us first note that for / G H SlyS2 we have 

\\-l- j\\H s i-- m i- rn '2' s 2- rn 3 S HiJ-si— mi-m2,S2-TO3-' /||l 2 

— I I Si - Wll— TO2I-S2 - J»3-* Si, — S2 I \L 2 ^L 2 I |-'-'-Si,S2 J I |i 2 

<cimi H 3i, S2 , 

if we show that n si _ mi _ m2jS2 _ m3 TIl_ Sl _ S2 is bounded in L 2 (IR n ). Assume (1). If we 
apply Theorems 12.11 and 12.51 (and also the argument of Remarks 12.21 and I2.6J) to the 

-SI,— S2 



composition (D) S2 m3 T(D) S2 , we can conclude that S = Il Sl _ mi _ m2iS2 _ m3 TIl_ 



is an operator of the form (13. ip with amplitude c(x, y, £) satisfying 

(3.6) \d^djc(x } y^)\ < C af37 (xr- m *(y) m2 - s ^\ 

for all a/3, 7, and for all x,y,£ G M™. Here we have used the fact that conditions 
(13. 5p guarantee conditions on the phase function in Theorems 12.11 and 12.51 (similar to 
the proof of Theorem 13.31) . Now, in the case (1), by Theorem 13.31 we obtain that S is 
bounded in L 2 (R"). 

Now assume (2). By Remark [2T71 we get that the amplitude of the operator S above 
is again of the form (13. ip and it satisfies estimate 

\d*d^c(x,y,0\ < C aPl (xr- m2 - H (yr- s \ 

again implying the Z/ 2 (R n )-boundedness of S by Theorem 3.1 from |RS4] . 

□ 
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4. Global smoothing for hyperbolic equations 

In this section we will consider the following Cauchy problem 

/ 4 1 n f id t u - a(D x )u = 0, 

^ ' ' y u(0,x) = u (x). 

The estimate of Theorem 14. II of this section can be also viewed as a quantitative way 
to express the fact that space decay of Cauchy data u translates into the time decay 
of the solution u. 

Let us assume that operator a(D x ) has a real valued smooth symbol a = a(£) e 
C°°(R n ) satisfying 

Vo(0 ^ for all f G E n . 

This condition makes the equation "dispersive", ensuring that the singularities will 
escape to infinity as t — > ±oo. The dispersiveness is in general necessary to have the 
smoothing properties for evolution equations, although some smoothing still takes 
place without this assumption (see |RS6] ) . 
Assume also that 

«(0 = oi (0 + o (0 for large f , 

where a\ = Oi(£) € C 00 (R™ \ 0) is a positively homogeneous function of order 1 (i.e. 
a x (A£) = Aai(£) for all A > and £ ^ 0) satisfying Va x (£) 7^ for all £ =£ 0, and 
a o = oo(0 G 5 10 is a symbol of order zero, i.e. a smooth function such that for every 
multi-index a there exists a constant C a > such that 

\d a a (0\ < C Q (0~ H foralieel". 

We note that from these assumptions it follows that a itself is a symbol of order one. 
In the case n — 1, we can easily see that the following estimate for the solution 
M = u(t, x) of (14. lp is true for s > 1/2: 

(4.2) H(a;)~ s «||L 2 (iR t xiRs) < C||mo||l 2 (rs)- 

Note that if a(£) itself is homogeneous, then w(t, x) is just a translation of m , so (14. 2 j) 
trivially holds. In |RS2j we have presented a result (which in turn follows from more 
general results proved in [RS6j ) that implies that estimate (14. 2p is also true for all 
n > 1. Analysis of previous sections allows us to apply this to obtain further global 
properties of solutions to (14.11) . 

Theorem 4.1. Let u(t,x) be the solution of the Cauchy problem (14.11) . Then for 
s > 1/2 and for all k = 0, 1, 2, . . . we have estimates 

(4.3) ||**(a:)~*(a;)~*«||ia(H*xK!f) < Ck\\(x) k u \\ L 2 {R n), 
where constant C^ is independent of Cauchy data Uq. 

Proof. For k = formula f!4.3[) is the same as (14. 2p . Let now k > 1. We want to show 
that estimate 

(4-4) \\t k (x)~ Sk u\\ L 2 {RtxR n) < C fc ||(x) fc M ||L2(K S ), 

holds for Sk > k + |. 
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The solution u = u(t, x) to (14. ip is given by 

u(t,x) = T t u (x)= f f e^ x -^ +ta ^u (y)dydC 



n ./ran 



Let a(X,D) be a pseudo-differential operator with symbol a(x, £) = x ■ Va(£). Let 
us denote 

r(£, x, := a(x + £Va(0, = (x + tVa(O) • Va(f). 
Taking 0(x, £) = x ■ £ + ta(^), it is straightforward to see that function \1/ defined in 
Theorem 12.81 is 

tf (77, £, ») = ^ • (£ - 77) + t(a(0 - 0(77)) + (77 - • (x + tVa(O). 

We can observe that V v e l9 ^' x ^\ v= ^ = 0. We can also observe that the corresponding 
function ty(x,y,£) in Theorem 12.51 is identically zero. Since a is affine in x, it follows 
from Theorems 12.51 and 12.81 that 

r(t,X,D)oT t = T t oa(X,D). 

Now, taking a(X, D)uq as the Cauchy data in (14.11) and using (14.21) with s > 1/2, we 
get 

1 1 (x)~ s T t o cr(X, D)u \ |L2(R t xRs) < C\ \<r{X, D)u \ | L 2 (B n). 
This means that we have the estimate 

\\(x)~ a (x + tVa(D)) ■ Va(D)u|| L 2 (RtxM n ) < C\\x ■ Va(D)«o||L3(R») 

for s > 1/2. Since Va(D) is a vector of pseudo-differential operators of order zero, 
we have the estimate 

\\x ■ Va(D)M ||L2(M S ) < C , ||(x)m ||l2(r S ) 

using the L 2 (M n )-boundedness of the operator x ■ Va(D)(x)~ by Theorem 13.51 Ap- 
plying the same argument k times, we get estimates 

(4.5) \\{x)~ s [{x + tVa{D)) • Va(D)] / M|| L 2 (RtXffi >n ) < C\\{x) 1 u \\l^(^), 

for all I = 1, . . . , k and s > 1/2. Now, it follows that 



\t k {x)~ Sk u\\ L 2 



(RtxKg) 
2k 



(4.6) =\\(x)~ Sk \Va(D)\- 2k l(x + tVa(D))-Va(D)-x-Va(D)} k u\\ L 2 (RtxR r i) 

< C\\(xy Sk [(x + tVa(D)) ■ Va(D) - x ■ Va(L>)] fe u|U 2(RtxRS) , 

where we have eliminated | Va(D) \~ 2k from the estimate using the L 2 (IR n )-boundedness 
of the operator (x)~ Sk \'Va(D)\~ 2k (x) Sk by Theorem 13.51 Now, we can estimate 

2 (M t xR™) 



(a;)-* [( x + tVa(D)) ■ Va(D) - x ■ Va(D)f u\\ L 2 {Rt> 



(4.7) 



k 

<Y, C kW( x )~ Sk i x ■ Va(D)] k - 1 \(x + tVa(D)) ■ Va(D)f u\\ LHRtXm) 
1=0 

k 
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where once again we used the weighted L 2 (IR n )-boundedness of \Va(D)\ k ~ l to elimi- 
nate it from the last inequality. We also note that operators in (14. 7p are of the same 
form, or are scalar, and they satisfy weighted L 2 (R n )-boundedness, so they have good 
commutator properties as well. Since 

s k -k + l>k + --k + l>-, 
we can use ( 14.51) . Combining it with f!4.6[) and (14.71) . we conclude that 

\\t k (x)~ Sk u\\ L 2 {RtxR n) < C||(x) fc U ||L2 (R n ) , 

completing the proof of (14. 3 p . □ 

We remark that estimate (14.31) with k = 1 is sharp in the sense that in general it 
does not hold for s = 1/2. For example, its failure can be easily seen for n = 1 and 

We have the following version of Theorem l4.1l if a(£) itself is positively homogeneous 
of order one. In this case a(£) has a singularity at the origin which causes a small 
additional complication when we eliminate operators of the form \Wa(D)\~ 2k from 
weighted estimates. In this case we have the following corollary from the proof of 
Theorem 14.11 

Corollary 4.2. If a(£) is positively homogeneous of order one with Va(£) 7^ for all 
£ 7^ 0, estimate (14.31) of Theorem\4-l\ still holds for < k < v ^- 



To justify it, we note first that operator |Va(D)| has symbol which is positively 
homogeneous of order zero. Consequently, for any I, operators (x) ~ Sk \'Va(D)\ l (x) Sk 
is L 2 (R n )-bounded for — ^ < s^ < | (see Kurtz and Wheeden |KW] . for example). 
Since we also need s k > k + |, we see that the conclusion of Theorem 14.11 remains 



valid for k < r± ^. 
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